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In order to derive the reciprocity relations, Onsager formulated a relation between thermal equi- 
librium fluctuations and relaxation widely known as regression hypothesis. It is shown in the present 
work how such relation can be extended to finite and isolated classical systems. This extension is 
derived from the fluctuation-dissipation theorem for the microcanonical ensemble. The results are 
exemplified with a nonintegrable system in order to motivate possible applications to dynamical 
systems and statistical mechanics of finite systems. 
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I. INTRODUCTION 

In his seminal work on the reciprocity relations [l[, On- 
sager formulated a hypothesis about the decay of ther- 
mal equilibrium fluctuations that was essential in his 
derivation. From this hypothesis, Onsager was able to 
relate the relaxation of nonequilibrium macroscopic ob- 
servables, obtained from phenomenological equations, to 
the decay rate of equilibrium fluctuations of those quan- 
tities. Such a relation between those two apparently dif- 
ferent phenomena is however quite expected from the 
point of view of the fluctuation-dissipation theorem de- 
rived [1] 20 years after Onsager's work. Indeed, since 
the derivation of the reciprocity relations from linear 
response theory, Onsager's regression hypothesis is un- 
derstood as a different but equivalent statement of the 
fluctuation-dissipation theorem. 

On the other hand, Onsager's reciprocity relations 
have already been extended to far from equilibrium con- 
ditions in different contexts [3] where nonlinear ef- 
fects are taken into account. Those extensions are based 
on fluctuation theorems which quite often are derived 
from the so-called "chaotic hypothesis" Q . Hence the 
fluctuation theorems play the role of linear response the- 
ory beyond near equilibrium conditions and have been 
also applied to understand the nonequilibrium behavior 
of systems far from the thermodynamic limit Q . 

Onsager's original work and its extensions mentioned 
above always focus on systems in contact with reservoirs. 
Here however we intend to study the near equilibrium 
behavior of isolated systems when they are finite and 
standard linear response theory cannot be applied. The 
relaxation to equilibrium of finite and almost isolated 
quantum systems has been studied experimentally us- 
ing cold atoms and has motivated several theoretical 
works which try to understand it and discuss the con- 
troversies related to the role played by the nonintegra- 
bility in this process (sl-fioj. For classical systems under 
such constraints, we show here that a simple extension 



of the fluctuation-dissipation theorem [ll|, li3 shows very 
clearly the importance of the dynamics in the relaxation 
to equilibrium within a description that is essentially 
Onsager's regression hypothesis extended to this new 
situation. An important extension of the fluctuation- 
dissipation theorem for the microcanonical ensemble was 
first derived by Nielsen [ll|. It deals with thermody- 
namic responses to thermal disturbances as for example 
heat pulses. In Ref . [l^ , only the response to mechanical 
disturbances, i.e., those which can be described as addi- 
tional terms in the Hamiltonian, are considered. Since 
they treat different aspects of the same subject, these 
works complement each other. 



II. DERIVATION 

We start presenting the usual regression hypothesis ex- 
pressed in mathematical terms. Let us consider a system 
described by the following Hamiltonian 



f)J-f 

H{Xo + dX) = H{\o) + dX — 

oX 



(1) 



A=A„ 



When the value of the parameter A is suddenly switched 
from Ai = Ao + dX to Ao at t = to-, the nonequilibrium 
av erag e value B of an observable B of the system evolves 
as [l^ 

Bit^to)^{B),^--^C,At-to), (2) 

for t > to, where {■)\^ denotes the equilibrium average 
value when A = Aq, fc^ is the Boltzmann constant, T is 
the temperature and 



CxAt) - {SA{0)SB{t))x, 



(3) 
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is the correlation function with 5X{t) — X{t) — {X)xi 
and A = {dH/dX)\x=x^. 

Equation 1^ is Onsager's regression hypothesis ex- 
pressed mathematically. It states that the relaxation of 
B to the equilibrium value is possible as long as the cor- 
relation function Cx^ (t) decays to zero. On the other 
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hand, if the relaxation process can be described by some 
sort of phenomenological equation, one obtains the decay 
of CAi(i) from ([2]). Since Eq. ^ is derived in the usual 
context of the canonical ensemble, the decay of the cor- 
relation function is interpreted as a consequence of the 
heat bath influence. 

We will now derive ([2|) in a different context, namely, 
when the system is isolated and finite, i.e. there is no 
heat bath and the number of degrees of freedom is such 
that the system is not in the thermodynamic limit. We 
consider first the system in equilibrium under H(Xo+dX). 
When t is equal to to, Ai is suddenly switched to Xo- The 
system then relaxes to a new equilibrium state. Assum- 
ing that the system was not far from the final equilib- 
rium state, a linear response calculation describes the 
relaxation process. Although the situation requires the 
microcanonical ensemble, the fundamental equations of 
linear response theory do not rely on any particular en- 
semble For t > to the Hamiltonian is H{Xo) and 
the situation can be stated as follows: a system, whose 
Hamiltonian is initially given by ([T]), is in equilibrium 
when at t — to the generalized force dX is suddenly re- 
moved and H{Xo + dX) —5- H(Xo). Therefore, from linear 
response theory, one obtains [lj| 

B{t-to) - {B)x, =dX f ds4)BA{E,t- s)Q{s~to), 



(4) 

where 8(2:) is the step function and {B)\^ is the following 
microcanonical average over the phase space points (q, p) 



with the distribution function p^i given by 
,5[i?-i/(q,p,Ai)] 



PAi(q,p) = 



ZxAE) 



where Zx^ [E) = / dq dp (5 [i? — i? (q, p, Ai 
The response function (pBAiE^t — s) is [1 



{E,t 



(5) 



(6) 



(7) 



where {•, •} is the Poisson bracket, 5X{t) = 
X(q(t),p(t)) - (X)ai, (q(t),p(t)) is the solution of 
Hamilton's equations of motion for N degrees of freedom, 
B is any observable and A = (dH/dX)\\=\„. 

The fluctuation-dissipation theorem [l2| yields then 



ZU! 



(8) 



where 



Xba{z,uj) 



^^g-(.^r+ij.) [Zx{E)4>ba{E,t)] 



Fba{z,u}) = 

— dr dEe-^^^^+^^-'^[Zx{E)CBA{E,T)], 

(10) 

and Cba{E, r), with t — t — s, is the following correlation 
function 



CBAiE.t^s) = {SAis)SB{t))x, 



(11) 



which differs from (jSj because of the microcanonical av- 
erage. 

Thus, Eq.® leads to 

^BAiE,r)^-^^^[ZxiE)CBAiE,r)]. (12) 
The integral in (|4]) can be written as 



dscpBAit - s)Q{s - to) ^ / drcpBAir), (13) 

CXD Jo 



and from (|12l) and (IT^ one obtains 



dT(f)BA{T) 



ZxdE 



[ZxCba{t)] 



(14) 



The dependence on E was omitted for convenience. 

Therefore, the relaxation of B to the new equilibrium 
state is described by the following expression for t > to 

B{t - to) = (B)a„ - dX^^ [Zx,CBA{t - to)] (15) 
Zai dE 



where 



1 d 



(9) 



{B)x^ = {B)x, + dX—— [Zx.CBAm , (16) 

since lim7._j.00 Cba{t) = is assumed. 

Analogously to Eq. ([5]), the relaxation to equilibrium 
of B is ruled by a correlation function related to it and 
the relaxation rate is given in terms of the decay rate 
of equilibrium fluctuations. Thus, the physical contents 
of Eq. (fT5|) allow us to interpret it as an extension of 
Onsager's regression hypothesis to the context of finite 
and isolated classical systems. 

Despite the analogy with Ea.([T5|) shows that 

the relaxation rate of B [as well as the decay rate of 
CBA{E,t)\ is given only by the statistical properties of 
the dynamics produced by H{Xo). Since the system is 
isolated, there is no influence of external thermal fluctu- 
ations on the decay of CBA(E,t) or B as in ([2]). There- 
fore, one might ask for which kind of dynamics relax- 
ation occurs. For integrable systems, whose motion is 
quasiperiodic, CBA{E,t) would decay only for N ^ 00. 
Nevertheless, it is well known that CBA{E,t) decays for 
chaotic systems (l5l - [l7| . For nonintegrable systems, the 
complete spectrum of behaviors, from quasi-periodic to 
chaotic, could be approximately obtained. 
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FIG. 1. Correlation function CBA{E\t) for B = A = 
{it + 92)/4 and E' = 2.5 The solid black lines were obtained 
numerically for 4x W initial conditions. The dashed red lines 
are the fitting of y4'e"°'* cos (oj't). (a) Ai = 0.1, A' = 21.5, 
a' = 0.20 and ui' = 1.7. (b) Ai = 0.12, M = 14.5, a' = 0.17 
and u)' — 1.8. 



In 1971, van Kampen made severe criticism of linear 
response theory ^18J which here, in the context of iso- 
lated and finite systems, seems to be even harder to an- 
swer. However, some of the arguments in the literature 
[13, H^, US] supporting standard linear response theory 
come from dynamical systems theory and are well suited 
for the present discussion. First, it is indeed possible, 
as mentioned above, that a finite and isolated system 
shows correlation functions which decay with time. In 
particular, if it has a statistical property called mixing 
[Tsl [l9| , it is possible to prove that such decay necessar- 
ily happens. Second, mixing is also responsible to ensure 
that an arbitrary smooth distribution in phase space will 
approach the microcanonical one for long times. Thus, 
although trajectories are very sensitive to small perturba- 
tions, the time evolution of distributions is rather stable. 
This would justify the linearization procedure leading to 
Eq. (|4]) at least for a class of systems. There should be 
also a constraint on time scales since the response func- 
tion (jl2[) would be ill defined for times much shorter than 
the inverse of the decay rate of correlations. 




FIG. 2. Relaxation of B for B = (Qi -1-92)74 and initial energy 
E = 5.0. The solid black lines were obtained numerically 
for 4 X 10^ initial conditions. The dashed red lines are the 



predictions of Eq. (|15p using Cba obtained numerically and 
the scaling property, (a) Ai = 0.1 switched to Ao = 0.12 and 
the value of (-B)ai = 7.54 obtained numerically for 4 x 10^ 
initial conditions with E — 5.0 and Ai — 0.1. (b) Ai — 0.12 
switched to Ao = 0.1 and the value of {B)x-^ — 6.55 obtained 
numerically for 4 x 10^ initial conditions with E = 5.0 and 
Ai = 0.12 



III. EXAMPLE 

In order to motivate possible applications of Eq. 
to both low-dimensional dynamical systems and statis- 
tical mechanics of finite systems, we will consider the 
relaxation process in a system described by the following 
Hamiltonian 1211 



H = 



9 

P2 



(17) 



which is integrable only for A = 1. It can be verified 
through the Poincare surface of sections that the motion 
generated by p?)) gets less and less regular as A decreases 
from the value of unity. If on one hand the relaxation pro- 
cess is very well defined for low-dimensional hyperbolic 
systems, on the other hand most of the realistic models 
used in statistical mechanics of classical systems are non- 
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integrable. For a nonintegrable systems with as few de- 
grees of freedom as p?)) , relaxation may or may not hap- 
pen and one has to find numerically the range of param- 
eters where it occurs. In our application of ([T5|) to this 
case, we have chosen A = {dH / dX)\\=\^ ~ [qf + (7|)/4 
and B = A. There is a small range of values of A 
(0-9 < A < 0.13) where Cba{t) can be fitted by the 
expression ^e"""^ cos (ojt) (see FigH]). It is not our aim 
here to find the exact functional form of the correlation 
function. Instead we want to find out whether a simple 
description of it (as the one just written above) is enough 
to describe the relaxation process approximately. 

The dynamics given by ([T7| is scalable with energy, 
i.e. the Hamilton equations remain invariant under a 
transformation {q[ 2,p'i 21^') ~^ (91,2,^1.2, i) given by the 
equations 



2hi=f£y^' pi^=(ex" 

9b \E') >i.2 UV 'V \e) ' 

(18) 

This property of (HH) yields A/ A' ^ {E / E'f since 
{qf + q|)2 oc {E/E'f, a/a' = {EjE'^l^ and ^juj' = 
{ElE'fl^. One can also verify that ZxiE) oc E^l"^. 
Therefore, the scaling with energy allows us to perform 
the derivative in dTSl) taking a certain value E' as the ref- 
erence and obtaining CBA{E,t) from the parameters of 
Cba{E' ,t)^ keeping of course the same functional form 
and the same value of A. In summary, the system de- 
scribed by p7|) was chosen as an example because it al- 
lows a simple illustration of (jlSp applied to extremely 
few degrees of freedom. For models with Lennard- Jones 
potentials, for example, the dependence of correlation 
functions with energy is much more complicated and is 
accessible only numerically. Besides, it has been shown 
recently [l^ that a finite collection of (|T7l) can act as an 
environment that induces relaxation on a simple degree 
of freedom. 



IV. DISCUSSION AND CONCLUSIONS 

All the numerical results were obtained from the inte- 
gration of the equations of motion of (1171) using a fourth- 
order sympletic integrator (23j . In Fig [1] is shown the nu- 
merical results for Cba{E' ,t). Although it is clear that 
the fitting is not excellent, one obtains afterwards a good 
agreement between numerical and analytical results for 
the relaxation of B. The analytical results should indeed 
be called semi-analytical since both {B)x-^ and Cba were 
obtained numerically. In Fig HI B[t — to) is the result of 
the average over several B{t — to) obtained from the time 
evolution of initial conditions distributed over an energy 
surface with E = 5.0. In Figl2Ka) a sudden switching of 
Ai = 0.1 to Ao = 0.12 at to leads to the relaxation of 
_B to a new equilibrium value. Although amplitude and 
frequency of oscillations are not correctly described, the 
relaxation time and the value of {B)\^ are well predicted. 
In Figl2Kb), the comparison between numerical and ana- 
lytical results for B is shown for Ai = 0.12 switched to 
Ao = 0.1. As before, amplitude and frequency of oscilla- 
tions are roughly described and the value predicted for 
(_B)a„ is not as good as in Figl^Ja). The relaxation time 
however is still in good agreement with the numerical 
results. 

In conclusion, we have derived an extension of On- 
sager's regression hypothesis from linear response theory 
when the system of interest is isolated and finite. Al- 
though thermal fluctuations induced by a heat bath are 
absent in this context, the expression obtained relates, as 
usual, relaxation to equilibrium fluctuations. Hence the 
new feature is that the decay of correlations is given only 
by the instrinsic dynamics of the system. The relaxation 
of a nonintegrable systems with two degress of freedom 
has illustrated that. The outlook is to extend the present 
result to the quantum regime. 
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